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Magnetic and charge noise have common microscopic origin in solid state devices, as described by 
a universal electron trap model. In spite of this common origin, magnetic (spin) and charge noise 
spectral densities display remarkably different behaviours when many-particle correlations are taken 
into account, leading to the emergence of the Kondo effect. We derive exact frequency sum rules for 
trap noise, and perform numerical renormalization-group calculations to show that while spin noise is 
a universal function of the Kondo temperature, charge noise remains well described by single-particle 
theory even when the trap is deep in the Kondo regime. We obtain simple analytical expressions for 
charge and spin noise that accounts for Kondo screening in all frequency and temperature regimes, 
enabling the study of the impact of disorder and the emergence of magnetic 1// noise from Kondo 
traps. We conclude that the difference between charge and spin noise survives even in the presence 
of disorder, showing that noise can be more manageable in devices that are sensitive to magnetic 
(rather than charge) fluctuations and that the signature of the Kondo effect can be observed in spin 
noise spectroscopy experiments. 

PACS numbers: 72.70.+m, 75.20.Hr 


I. INTRODUCTION 

The tunneling of conduction electrons into local 
charge traps is a prevalent phenomena in solid state 
physics. Traps can be realized by artificial structures 
such as quantum dots,® or by natural “unwanted” de¬ 
fects such as dangling bond:P and bound states in 
metal/oxide interfaces.® It has long been recognized that 
trap fluctuation causes charge noise in electronic devices, 
with the signature of individual traps being observed 
with a L oren tzian l// 2 noise spectral density in small 
structures® and an ensemble of them causing 1// noise 
in large structures. 6 Here we address the fundamental 
question of how the electron spin alters trap noise. 

One of the greatest developments of interacting elec¬ 
tron physics was the discovery that a local trap inter¬ 
acting with a Fermi sea gives rise to the Kondo effect, 
the formation of a many-body singlet with conduction 
electron spins screening out the local trap spin.® The sig¬ 
natures of the Kondo effect in transport phenomena are 
well studied, but key issues related to dynamics have only 
been addressed recently with the emergence of modern 
Numerical Renormalization Group (NRG) algorithms.® 
It is particularly interesting to And out whether trap 
noise will impact devices that are sensitive to magnetic 
fluctuations a s oppo sed to charge , e.g. spin-based or spin- 
tronic devices,®® 1 in the same way that it affects con¬ 
ventional charge-based devices. Recent measurements of 
intrinsic magnetic flux noise in superconducting quantum 
interference devices do indeed confirm that trap spin fluc¬ 
tuation is the dominant source of noiseP®^ Moreover, 
novel developments in spin noise spectroscopy 15 open 
several possibilities for the detection of correlated spin 
fluctuations in quantum dot systems. 


Given these interesting prospects, the question that we 
address here is the qualitative difference between pure 
charge/spin noise of a “Kondo trap” interacting with a 
Fermi sea, which we define as a local charge trap in the 
Kondo regime. 

The interplay of Kondo physics and noise has been 
mostly explored in the context of transport through 
quantum-dot systems, with the Kondo trap right inside 
the transport path. In this case trap charge and spin 
fluctuation are intertwined in a non-trivial way. Calcula¬ 
tions of the shot noise and current noise in different set¬ 
ups such as singl(h-2H22J anc j double quantum dotJ^SEU j n 
the Kondo regime have been reported. Much less stud¬ 
ied is the role of the Kondo state in spin noise. The 
case of spin-current noise was considered in Refs. 1241 and 
1251 and qualitative differences between spin-current and 
charge-current noise were found to exist. 

In this article, we show that focusing on pure 
spin/charge trap noise (i.e., finite frequency trap occu¬ 
pation noise) allows for a different perspective on the 
problem of Kondo trap dynamics: it enables a clear sep¬ 
aration between the contributions of single-particle exci¬ 
tations and the many-particle processes connected with 
the formation of the Kondo singlet state. Moreover, con¬ 
sidering pure spin (charge) trap noise is important for 
describing transport experiments with traps outside the 
transport channel. In this case, trap fluctuations produce 
bias magnetic (electric) noise that in turn may dominate 
the spin-current (charge-current) noise. 

Our article is organized as follows. In Section [XT] we 
outline our model for pure spin/charge trap noise, and 
establish its connection to the usual spin/charge suscep¬ 
tibilities. We demonstrate six exact results: four sum 
rules and two Shiba relations. In Section nn we de- 
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scribe our Hartree-Fock (HF) or mean-field approxima¬ 
tion, that mainly accounts for single-particle processes. 
In Section [IV] we present our non-perturbative NRG cal¬ 
culations, which account for single-particle and many- 
particle processes on the same footing. The NRG results 
show that finite-frequency spin/charge noise have quite 
distinct behaviors and are dominated by completely dif¬ 
ferent processes. In Section [V] we use NRG and the sum 
rules to obtain an analytic approximation to spin noise 
in the Kondo regime, and in Section |VI| we use this an¬ 
alytic approximation to study the interplay between dis¬ 
order and Kondo correlations in an ensemble of Kondo 
traps. We show that, in the presence of disorder, the spin 
noise displays a temperature-dependent 1 // noise that is 
qualitatively distinct from the temperature-independent 
charge 1 // noise. Finally, Section VII presents our con¬ 
cluding remarks, with a discussion of the impact of our 
results in the effort to detect Kondo correlations in spin 
noise spectroscopy experiments, and our prediction of 
qualitatively different 1 // noise impacting spin-based 
and charge-based devices. 


II. CHARGE TRAP MODEL AND EXACT SUM 
RULES 

Our starting point is the Anderson modePS for a 
trapping-center interacting with a Fermi sea, 


(n-f+riq), respectively, and (•) denoting the thermal equi¬ 
librium average. 

We write an exact expression for the spin and charge 
noise by performing a spectral decomposition of Eq. © 
in the basis of energy eigenstates: 


^ 


i—Em/T 



5(u)-E nm )-(6i) 2 S(uj) 


( 4 ) 

where Z is the partition function, jm) are (many-body) 
eigenstates of the Hamiltonian (jT]) with energy E m 
( E nm = E n — E m ) and (n\Oi\m) are the many-body 
matrix elements of the local operator Cfo For simplic¬ 
ity, we set h = ks = 1- Note that Eq. © implies that 
Si(oj,T ) > 0 and Si(—iv,T ) = e~ ul / T T) as required 
by our assumption of thermal equilibrium. 

The noise spectra is closely related to the dynami¬ 
cal susceptibility associated with the operator Cfo We 
shall explore this connection in order to derive the exact 
frequency sum rules and Shiba relations^ for Si(u,T). 
These relationships will be used in Sec. [V] to obtain ana¬ 
lytical approximations for the noise spectra. 

Assuming that an external field Fi(t) couples to 
Oi through "Hext = — OiF^t), the linear response 

of Oi to Fi will be (O l (t)) Fl t 0 - ( Oi) F=0 = 
2n J duje~ lwt Xi(^,T)Fi(uj), where Xi(w,T) is the dynam¬ 
ical susceptibility given by2S 


H = Hh&nd + 'Hhyb + %trap, 

( 1 ) 

^band — ^ ^ ^ka^kcr 5 

( 2 a) 

^hyb — ^ ^ Vdk (^'ka^ J(T ^a^ka^ 5 

k . <7 

( 2 b) 

"Htrap = (n t + nj.) + Un t n±. 

( 2 c) 


In the above, c' k<T ( Cka ) is a creation (destruction) op¬ 
erator for a conduction electron with wavevector k and 
spin a =t, rika = c^Cka counts the number of band 
electrons in state k, a with energy eka- Similarly, the 
operators dj. and d a create and destroy a trap electron 
with spin cr, respectively, with n a = d\d a being the num¬ 
ber operator for electrons with spin cr occupying the trap 
state with energy e^. Finally, U is the Coulomb repulsion 
energy for the trap, with ed + U the energy required to 
add a second electron to an trap site that already con¬ 
tains one electron. 

Our goal is to calculate the trap spin S s (co,T) and 
charge S c (uj : T) noise spectral densities, defined by: 


Si= s ,c{^,T) = — 


dt e iut (50^)50^0)) , (3) 


W(^T) = - 


dt e l ‘ 


di(t),di( o) 


( 5 ) 


Performing a spectral decomposition of Eq. ([5]) and 
comparing to Eq. © leads to the following Lehmann 
representation: 


1 r°° rl , 

XMT) = / -K(-u/,T) - Si(u)',T )], 

7_ 00 w - Ld' + IT) 

( 6 ) 

with 77 —> 0 + . Separating the susceptibility into real 
and imaginary parts, \i = Xi + *x". using Si(-u,T) = 
e-“/ T Si( ojjT), and taking the imaginary part of Eq. © 
leads to 


Xi{u,T) = 1 ^ Si(u>,T), (7) 


which is known as the fluctuation-dissipation theorem. 
Moreover, taking the real part of Eq. © yields 


Xi(u,T) = V r -p— (l-e- w '/ T ) Si(u/,T), 
2 tt 7-00 u' -u \ ) 


( 8 ) 


which is the Kramers-Kronig causality relation. 

We now derive the frequency sum rules. The first one 
is obtained by direct integration of Eq. Q over all fre¬ 
quencies: 


where SOi(t) = Oi(t) — (Oi) with trap spin and charge 
operators given by O s = S z = (n-(- — 714.)/2 and O c = 


Si(oj,T ) dui 


( 6 2 ) - ( 4 ) 2 . 


( 9 ) 
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FIG. 1: (color online) Charge noise as a function of frequency 
for the trap in the symmetric case with = —U /2. NRG cal¬ 
culations are shown to be well approximated by a mean-field 
Hartree-Fock decomposition (HF) even when U/Y is large and 
the trap is deep in the Kondo regime. This shows that charge 
noise is well described by single-particle excitations. 


We call this the spin or the charge sum rule depending 
on whether i = s or i = c. Another sum rule is obtained 
by setting w = 0 in Eq. ([8]), and noting that Eq. ©> 
implies Xi(oJ = 0 ,T) = x'(w = 0 ,T): 


1 - e~ u ’'/ T 
27Tu/ 


Si(u,T)du' = X i(v = 0,T). (10) 


Accordingly, we call this the spin or charge susceptibility 
sum rule. Altogether Eqs. ([9]), (10) form a set of four 
exact sum rules that are valid at any temperature T. 

Finally, there are two additional exact relationships 
between noise and susceptibility, that apply only at T = 
0. These are the so called Shiba relations: 27 


Lim^ 0+ Ss ^ 0 'T_ °) = [ Xs ( w = 0 ,T = 0)f ,(lla) 


Lim, 


■ui— > 0 + ' 


87 r 2 w 
S c (uj,T = 0) 

27T 2 W 


= [x'c(w = 0, T = 0)] 2 (lib) 


They imply that Si(u), T) is Ohmic (linear in w) at T = 0, 
with a slope related to the static susceptibility Xi( w = 
0, T = 0). 


III. HARTREE-FOCK APPROXIMATION 


FIG. 2: (color online) Spin noise as a function of frequency 
for the trap in the symmetric case with td = —U/2. The 
NRG results agree with HF only at U = 0. As U increases 
NRG shows that the magnetic noise increases, developing a 
peak at u « Tk■ In contrast, the single particle contribu¬ 
tions described by HF decrease dramatically as U increases. 
This shows that magnetic noise is dominated by many-body 
processes. 


charge noise is 



( 12 ) 

and for the spin noise we get simply S^ F (ui,T) = 
|S , ^ IF (u;, T), i.e., in the HF approximation magnetic 
noise is simply j times the charge noise. In Eq. (12) 
/(e) = l/[exp ((e — ep)/T) + 1] is the Fermi function, 
and 


(e) 

Au.( e ) 


r / 7r 

(e- £d) 2 + T 2 ’ 

r / 7 t 

(e-e d -C/) 2 +T 2 ’ 


(13a) 

(13b) 


are HF local density of states for the trap with spin / and 
/, respectively. The energy scale T = 7r pV^ models the 
rate for escape of a trap electron into the Fermi sea, with 
p the energy density at the Fermi level, and Vdk = Vd&k- 
independe nt co uplin g bet ween trap and Fermi sea. Note 
that Eqs. (13a) and (13b) break the local spin symmetry 
by assuming the energy for the f and / trap states are 
td and ed + U, respectively. This result is well known 
to be incorrect, in that it misses Kondo physics, i.e. the 
screening of trap spin by the electron gas spins. 


As a first approximation we calculate the noise spectral 
densities using Hartree-Fock (HF) decomposition based 
on writing expectation values into products of spectral 
functionsPSJ The advantage of HF is that it becomes ex¬ 
act in the U = 0 non-interacting limit!^® The result for 


IV. NRG CALCULATIONS 

We shall compare the Hartree-Fock approach to non- 
perturbative NRG calculations of the noise spectra, that 
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FIG. 3: (Color online) Universal scaling for spin noise in the 
Kondo regime for td = —U/2. (a, b) NRG results for spin 
noise S s {lj). Note how all the curves collapse into a single 
scaling relation when the noise is written as a function of 
oj/Tk . For u> < Tk, the magnetic noise scales linearly with 
w (Ohmic noise), and for Tk < uj < U it decreases with 
an anomalous power of frequency oc l/[wlog 2 (w/Tk)]- For 
u>U, spin noise is cut-off oc 1/uj 2 . (c, d) NRG results for 
charge noise S c {uj) do not show universal Kondo scaling , and 
behaves just like the single particle approximation (HF) with 
noise peaked at u « Max{F, U} with smooth cut-off 1/u) 2 at 
u > U. 


take into account local spin symmetry and the formation 
of the Kondo singlet. The NRG algorithm calculates, 
within some well-controlled approximat ions, 9 the many- 
body spectrum for the Anderson modelP^D Conduction 
electrons are assumed to have a continuum spectrum, 
forming a metallic band with a half-bandwidth D. 

At zero temperature, the first term in Eq. © can be 
computed from the NRG spectral datcPSHEll down to ar¬ 
bitrarily small non-zero frequencies |w| > 0. The spec¬ 
tral weight at u> = 0 and the fulfillment of the sum-rules 
can be obtained by calculating the expectation values 
(Oi) and {Of) with NRG. Since we will be interested 
in the large frequency regime and our spectral functions 
obey well-defined sum rules, we have ch osen t o use the 
“Complete Fock Space” (CFS) approac H 32 * 331 to calcu¬ 
late S t (oj > 0) at zero temperature. As discussed in 
Appendix [Aj this choice has two important features: (i) 
the T = 0 spectral functions are sum-rule-conserving by 
construction and (ii) broadening artifacts in the high fre¬ 
quency regime, which can mask the correct power-law 
behavior, are minimized. . 


Figure [l] shows the calculated charge noise in the case 
ed = —U/2 for F = 10 ~ 4 D and several different U. Re¬ 
markably, HF remains a good approximation to charge 
noise even at large U. We interpret this result to be ev¬ 
idence that charge noise is dominated by single particle 
processes even when the trap is deep in the Kondo regime 
(U >T forT = 0). 

The situation is drastically different for magnetic noise 
as shown in Fig. [2] While NRG and HF agree with each 
other in the (7 = 0 limit (when HF is exact), as soon 
as U becomes non-zero the two methods show opposite 
results. As U increases, the single particle noise (HF) 
decreases, while the many-body noise (NRG) increases. 
The low-frequency NRG results can be better visualized 
in Fig. [3j We find (Fig. |3]-a,c) that the magnetic noise 
for a single trap is Ohmic at low frequencies, with a 
peak at w « Tk where Tk is the Kondo temperature. 
The magnetic noise spectral densities all collapse in the 
same universal curve and scale as an anomalous power 
law oc Tk/\lo log 2 (i u}/Tk)\ in the Tk -C w -C U frequency 
range (Fig.[3]-b), consistent with previous results for the 
dynamical spin susceptibility 34 38 and the spin-current 
nois^S in the Kondo regime. 


V. ANALYTICAL APPROXIMATION FOR 
SPIN NOISE IN THE KONDO REGIME 


While the HF approximation [Eq. ©] failed to de¬ 
scribe spin noise, it was shown to give a good description 
of charge noise at T = 0 [Fig. [I]. In Appendix [B| we 
show that HF actually provides a good approximation for 
charge noise at T > 0, in the sense that it approximately 
satisfies the sum rules and Shiba relations demonstrated 
in Section [Hj The goal of the current section is to use our 
NRG calculations, sum rules and Shiba relations to ob¬ 
tain an analytical approximation for spin noise at T > 0 
in the Kondo regime. 

It is well knowrPthat NRG has difficulty in calculating 
spectral features at frequencies u <T. Here, we propose 
an alternate approach to evaluate the spin noise for a 
broader u>/T range. 


Motivated by the susceptibility sum rule Eq. (101 and 
the property S s {—u,T) = e~ u / T S s (w, T) we propose the 
following fit function 


S™(u,T) = 


2wy s (w = 0, T) 


1 — e~ ul / T 


+ iT 


(14) 


with the u> = 0 susceptibility given by a continuous fit to 
the NRG resullP^ 
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Xs(uj = 0, T) 


1 

8 ttT 


1 


W 

8ttTk 

0.68 


8ti-(T+v / 2T a -) ’ 

i _ toefi<g{y/av» 

log (T/T k ) I log 2 (T/Tk) 


for T < 0.23 T k , 
for 0.23 T k <T< 28.59 T K , 

for T > 28.59 T K , 


(15) 


where W = 0.4128 is the Wilson number. 

In Eq. p| T s = T s (u,T) is a fit function of fre¬ 
quency and temperature that will be determined by the 
exact sum rules and the Shiba relations. We recall that 
previous relaxational fits for T, assume no frequency 
dependence.^ Here we allow r,(w,T) to vary on fre¬ 
quency so that the logarithmic frequency decay discussed 


in Section IV is properly accounted for. 

For ^ (V Tk, the perturbative method of Suhl- 
Nagaok d 41 ! 42 ! yields the high temperature limit (the Ko- 
rringa law): 


1 


r s (w, T » T k ) « — 
47T 


T 


1 + 


37r 


log 2 ( 


T 

T k 


(ie: 


At T = 0 the Shiba relation (11a) applied to Eq. (14) 
impliepSl 


T s (w = 0,T = 0) = 


2 T k 


47t 2 x s ( 0,0) 7rW 


(17) 


where we used the NRG result x«(0, 0) = W/(8ttTk). 

In order to interpolate between Eqs. ( fl6| and (17) we 
propose the following expression: 


r a (w,r) = 


l 


T+ %T k 


47r 


1 + 3^ lQ g 


1 + {tt) + (sfe) 


C 18 ) 

where a is a fit parameter to be determined by the spin 
sum rule [Eq. Q]: 


Sum, (T) = 4 / dwSZ 11 (u, T). 


(19) 


This sum rule is most sensitive to a at T = 0, and 
we find that the optimal fit value is quite close to a = 
3, when Sum, (T = 0) = 0.9994. As an independent 
check, we evaluate the spin sum rule at T > 0 and the 
susceptibility sum rule at T > 0: 


Sum Xs (T) 


i 


Xs(0, T) J_ c 


1 — p-^/T 

dco— -S s Fit ( W ,r). 

Zttlu 


( 20 ) 

In all cases, we obtain agreement within 36%. A few 
examples are shown in Table [T] Moreover, we find that 
Eqs. (14) and (18) with a = 3 provide an excellent fit of 
our NRG results at T = 0, as shown in Fig. [4j 

Note that the choice of Eq. ( fl8| implies that 
TkSJ^^jT) is a universal function of lj/Tk and T/Tk, 



FIG. 4: (C olor onli ne) Comparison of the spin noise fit 

iSf rt (w) [Eqs. (141 and (181 with a = 3] (lines) with the NRG 
results (symbols) at T = 0. 


T/T k 

Sum, 

Sum Xe , 

0 

0.9994 

0.9518 

0.5 

0.9247 

0.9502 

1 

0.8245 

0.9503 

10 

0.6910 

0.9875 

100 

0.7777 

0.9979 


TABLE I: Sum rules [Eqs. (191 and ( |20|) [ applied to our ana¬ 
lytical fit of spin noise, Eqs. m and ( |18p with a = 3. For the 
spin sum rules we used analytical approximations for y s (0, T ) 
obtained by NRG [Eqs. (4.53) and (4.60) in Ref. 2S1- In all 
cases we find that the sum rules are satisfied within 30%. 


and that the presence of the temperature-dependent 
functions x,(0,T) and T,(aj,T) suggest that spin noise 
has a much stronger temperature dependence than 
charge noise. In particular, Eq. (14) fully accounts for 
the Kondo screening for T<Tk through y s (u; = 0,T). 


VI. SPIN NOISE IN THE PRESENCE OF 
DISORDER 

In the case of an ensemble of N Kondo traps, the noise 
will be affected by disorder. The usual model for trap dis¬ 
order (the one that gives rise to ubiquitous charge 1/f 
noise)P is to assume trap tunneling rate T = T 0 e _A , 
where A models the tunneling distance between trap 
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(a) Low temperature behavior. 


(b) High temperature behavior. 


FIG. 5: (Color online) Spin noise in the presence of trap disorder. The calculated noise for N traps was averaged according to 
the prescription T = Toe~ A , with A the tunneling distance between trap and the Fermi gas uniformly distributed in the interval 
[0, A m ax]. This gives rise to the broad distribution of Kondo temperatures shown in Eq. ( |’22[ l. The resulting noise, shown here 
for k = 10 and A max = 5 displays 1// behavior over a frequency range that decreases as the temperature increases. This is in 
contrast to the temperature independent charge 1 // noise described in the literature.® 


and Fermi sea. The model assumes A uniformly dis¬ 
tributed with density P'( A) = IV/A max for A € [0, A max ], 
and P'( A) = 0 for A outside this interval, resulting in 
P(r) = (iV/A max )/r and the corresponding 1// fre¬ 
quency dependence for trap charge noise. As we shall 
show, this same model applied to Kondo traps gives rise 
to a much broader distribution of Kondo temperatures 
that we denote P(T K ). 

For definiteness, we assume all Kondo traps have fixed 
ed and C/, with the disorder solely affecting the parameter 
r(A). The dependence of the Kondo temperature with A 
is given bj® 


since the logarithm in Eq. (22) makes P(Tr-) flatter than 
a ~ 1/Tk distribution, thereby making the integral fi¬ 
nite. We remark that our P(Tk) is appropriate to de¬ 
scribe highly disordered traps, such as traps randomly 
distributed at an insulator close to the metal/insulator 
interface. This situation is quite different from Kondo 
impurities in bulk alloys, whose P{Tk) is considerably 
less broad .E311] 


Applying this averaging prescription to our spin noise 
Eq. (14) yields 


Tk( A) = 


/r(AW V3e d (, d + V) ! 

\ e u r ( A ) 

<S») = r 

V 2 t r 

JT\ 


dT K P{T K )Sg{uj,T). (23) 


= Tjy ax e _ [^ +K ( eA_1 )]. 


( 21 ) 


Here k = — \/3 ed{(-d + U)/(UT 0 ) > 0 characterizes the 
type of trap. We shall assume k (A max + l)/2, a limit 
that is typically satisfied by Kondo traps with (/>!. 
The maximum and minimum Kondo temperatures of the 
distribution are given by T™ ax = Tk {A = 0) and T™ n = 
T k { A = A max ), respectively; for T K e [T™ n ,T£ ax ] the 
trap density becomes 


P(T 1 - ® 

I cLT k 


N 


l dx\ Tk re-log (ySfir) 


( 22 ) 


with P(T k ) = 0 for T k </L [T“ in ,T“ ax ], Note 
how P{Tk) is exponentially broader than P(T): we 
have Tjy a,X /T™ n « exp [«exp (A max )], in contrast to 
k'max/k'min = exp(A max ). In spite of this difference, 
the normalization condition f cITkP{Tk) ~ N still holds 


The results are shown in Fig. (5}a,b. At low tem¬ 
peratures (T < T£? ax ) the noise shows 1// behavior 
up to frequencies of the order of T)^ ax ; at larger fre¬ 
quencies, Kondo-enhanced exchange processes lead to a 
l/[/log 2 (/)] behavior. For higher temperatures (T > 
jimax) t ] le no j se saturates in the low frequency region, 
and the l/[/log 2 (/)] behavior gets washed out of the 
high frequency region. 

Interestingly, the frequency range with 1// be¬ 
havior gets reduced as the temperature increases. 
This shows that spin 1 // noise behavior is strongly 
temperature-dependent, in marked contrast to the usu¬ 
ally temperature-independent charge 1// noise. The ad¬ 
ditional temperature dependence implies that tempera¬ 
ture actually competes against disorder, converting the 
spin 1 // noise into a Lorentzian. 
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VII. CONCLUDING REMARKS 

In conclusion, we presented a theory of charge and spin 
noise of a Kondo trap interacting with a Fermi sea. We 
showed that trap spin noise is qualitatively different from 
charge noise, in that the former occurs due to many-body 
scattering processes, while the latter is mainly dominated 
by single-particle tunneling. This difference implies that 
spin noise has a stronger temperature dependence than 
charge noise, and that it is controllable by tuning Kondo 
temperature Tk rather than trap tunneling rate T. 

Kondo trap dynamics displays two quite distinct be¬ 
haviors depending on which property is probed. The 
experimental methods of charge 1 ^ and spin 15 noise spec¬ 
troscopy use optical absorption to detect noise via the 
fluctuation-dissipation theorem [optical absorption at 
frequency w is directly proportional to x"(w,T) and to 
noise as in Eq. 0 ]- Our results elucidate how Kondo 
correlations can be observed with these methods. Pure 
charge absorption does not enable the detection of the 
Kondo effect; in Ref. |T| the formation of the exciton state 
mixes charge and spin fluctuation, and this feature was 
critical in enabling their observation of the Kondo effect. 
For spin noise spectroscopy, universal scaling with Ohmic 
behavior at T <u<T K coupled with a l/[wlog 2 (lo/T k )\ 
tail for Tk can be taken as the signature of the 
Kondo effect, allowing the extension of this technique 
to probe Kondo correlations. However, in the presence 
of strong disorder over a range of Kondo temperatures 
Tk € [Tjj? in , Tjj? 8 *], we find that the Ohmic behavior is 
washed out, and the signature of Kondo correlations are 
visible only for to > T“ ax and T < 10T“ ax [See Fig. [HJa]. 

The qualitative difference between spin and charge 
noise survives even in the presence of disorder and high 
temperatures (namely Tk T]^ ax ). As the tempera¬ 
ture increases, the range of 1 // behavior for spin noise 
decreases, while the range of 1// charge noise remains es¬ 
sentially unaltered. The additional temperature depen¬ 
dence for spin noise implies that temperature actually 
competes against disorder, converting the spin noise 1 // 
behavior into a Lorentzian-like dependence. Given that 
1 // noise is notoriously difficult to control ^ we reach 
the conclusion that ubiquitous trap noise can be more 
manageable in spin or flux-based devices that are sensi¬ 
tive to magnetic fluctuations rather than charge. 
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FIG. 6: (Color online) NRG data for the spin noise S s (lo) 
at T = 0 for U = 40r calculated with DM-NRG and CFS 
procedures. Inset: A check of the spin sum rule given by Eq. 
|Al| for the different approaches shows that the CFS fulfils the 
spin sum rule apart from numerical integration errors. 


Appendix A: Details of the NRG calculations 

As we argue in the main text, our choice of the Com¬ 
plete Fock Space (CFS) procedure^ (or, equivalently, the 
full density matrix NRG method (FDM-NRG)P^at T = 0) 
in the NRG calculations presents some advantages for the 
calculations of the correlation functions listed in Eq. Q . 
To illustrate this point, we compared results obtained 
using CFS and the earlier “Density Matrix-NRG” (DM- 
NRG) methodP 

The main panel in Fig. [6] presents NRG data for the 
spin noise S a (w) using DM-NRG (open circles) and CFS 
(filled squares) for U = 40r and other parameters set 
as in Fig. [3| In both cases, the NRG calculations were 
performed using a discretization parameter A = 2.5 re¬ 
taining up to 1000 states at each NRG step, which en¬ 
sures convergence for the single-trap Anderson model. 
The spectral data was broadened using the usual loga¬ 
rithmic Gaussian functions (Eq. (74) in Ref. EJ) with a 
broadening parameter b = log(VA) ~ 0.46 (We have used 
z-averaging for some of the data presented, particularly 
the data presented in Fig. [I]). 

Clearly, DM-NRG subestimates the peak at ui = Tk in 
comparison with CFS. More importantly, it misses the 
transition from the S s (Tk oc l/[wlog 2 (co/Tk)\ 

to S s (qj U) oc w -2 behaviors, which is one of the 

important features distinguishing the spin noise from the 
charge noise. 

We have also compared both methods by checking the 
fulfillment of the spin sum rule in Eq. §■ In the absence 
of spin polarization (due to, e.g. magnetic fields or fer¬ 
romagnetic couplings), (S z )=0 and the spin sum rule is 










given by: 


S s {u J )duj = (S 2 z ). 


(Al) 


To this end, we performed a static NRG calculation of 
( S Z )(T —> 0) (open squares in the inset of Fig. |6j) and 
compared with a numerical integral of S s (ui). The agree¬ 
ment of the integrated CFS data (diamonds) is much 
better than the DM-NRG (circles), although the fulfill¬ 
ment of the spin sum rule is not perfect due to numerical 
errors coming from the integration procedure. 

Alternatively, the integral in Eq. (Al) can be done di¬ 
rectly in Eq. Q and it becomes a sum over the spec¬ 
tral weights | (n \ S z \ m) \ 2 provided that the set of many- 
body states {|to)} forms a complete set. In practice, this 
procedure can be done in the CFS scheme, as it retains 
matrix elements between “discarded” and “kept” NRG 
states, maki ng th e set of many-body states complete by 
construction. 32 — In this case, free of numerical integra¬ 
tion errors, the CFS data (filled squares) fulfills the spin 
sum rule down to machine precision, as shown in the inset 
of of Fig. § 


Appendix B: Validation of Hartree-Fock 
approximation for charge noise when T > 0 


The HF approximation Eq. (12) was shown to approx¬ 
imate charge noise at T = 0. Here we check its validity 
at T > 0, by direct evaluation of the Shiba relations and 
sum rules described in Section El 

The static (w = 0) charge susceptibility in the HF 
approximation is given by 


,HF 


( w = 0 ’ T ) = 


27t de F 

1 


de}2Ava(e) g ^ 

_ [ de ^ A(Tcr ^ 

87 tT J cosh 2 ( 5 =|K| 


At T = 0 we get 


(Bl) 


HF 


Xc 


( W = 0 ,T = 0) = — £>„(€*■). (B2) 


We start by checking the Shiba relation for charge 
noise, Eq. (lib). In the HF approximation we get 


reF+ui ^ 

Lim w _> 0+ / — Y A <7 „(e)A acT (e-uj) = ’V'A 2 CT (e F ), 

J<LF W 

a <7 

(B3) 

which according to the Shiba relation should be equal to 


2n 2 [ X ? F (0,0)] 2 =l 


YA aa (e F ) 


(B4) 


T/r 

e d /r 

U/Y 

Sum Xo 

0.1, 

1, 

10 

0 

0 

0.9997, 

0.9999, 

1.000 

0.1, 

1, 

10 

-2.5 

5 

0.9994, 

0.9999, 

1.000 

0.1, 

1, 

10 

-10 

20 

0.9915, 

0.9978, 

1.000 

0.1, 

1, 

10 

0 

5 

0.9996, 

0.9999, 

1.000 

0.1, 

1, 

10 

-10 

10 

0.9996, 

0.9999, 

1.000 


TABLE II: Charge susceptibility sum rule in the HF approx¬ 
imation [Eq. ( |B6| )]. The sum rule is seen to be satisfied 
(Sum Xc = 1) with high accuracy for several different param¬ 
eters. 


The relation is satisfied exactly at U = 0; however, as U 
increases Eq. (B3) becomes up to two times larger than 


Eq. (B4). This discrepancy can indeed be observed in 


the comparison with NRG, see the difference in slopes at 
u> = 0 in Fig. |T] Nevertheless, the discrepancy is not too 
large. 

The charge sum rule [Eq. ([£]) for * = c] in the HF 

flnnrnYimflt.irm rparls 


f dojS? F (u,T) = W deA a <j(e)[l — /(e)] 

x J didA aa (e - w)/(e - w) 


= - ( n °)]( n °) 
a 

= {6 2 c ) hf -(6 c ) 2 hf, (B5) 


where ( O c ) HF = f deA ac7 (e)f(e) and (6 2 ) HF is ob¬ 
tained by making t he a pproximation (n-j-n^) ~ (n^)(ni). 
The last line of Eq. ( |B5[ ) is expected to be a good approx¬ 
imation to the exact result, even in the Kondo regime, 
when charge fluctuations are strongly suppressed. 

Finally, we verify the charge susceptibility sum rule 
[Eq. ( |10[ ) for i = c] with explicit numerical calculations 
of the quantity 


Sum v = ———- / 

x ? F {o ,t)7_ c 


f — p -ui/T 

du— - S ff (co,T). 

Z7 TUJ 


(B6) 

As shown in Table [H] these values are very close to 1 for 
all tested parameters. 


In conclusion, the HF approximation for charge noise 
is consistent with the exact relations of Section El for all 
parameters checked, indicating that it provides a good 
analytical approximation for charge noise even for T > 0. 
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